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The phase-space paths introduced by Cheuk-Yin Wong in Phys. Rev. C25, 1460 (1982) and 
discussed recently in the literature can be used for calculation of evolution of the Wigner function 
to first order in the time increments only. The first-order solutions are helpful to determine the 
phase-space Green function in the framework of the phase-space path integral method. 
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^ : 

Quantum trajectories exist in the de Broglie - Bohm theory [l], [3, S| and appear in the framework of the deforma- 
tion quantization i4 i j, B , H, @ as the Weyl's symbols of the Heisenberg operators of canonical coordinates and 
O ^ momenta [lol . [ill . lUllJ • In the de Broghe - Bohm theory, the particle trajectories play an important role in the 
1^ interpretation of measurements, whereas within the deformation quantization framework quantum trajectories have 
C/^ , properties assigned to characteristic lines of first-order partial differential equations (PDE). The evolution equation 
— . ' for the Wigner function is the infinite-order PDE, nevertheless, it can be solved with the help of quantum character- 
' istics. The interest to particle trajectories is motivated by the fact that they remain striking but intuitive feature of 
transport models in heavy-ion collisions 



Already some time ago, C. Y. Wong 



18,|19|,|20|,|2H. 

Jj| proposed a new kind of quantum trajectories for solving the evolution 
^ _ problem of the Wigner function. Recently these ideas were discussed in Refs. [l6l[l7|. 

I ' We wish to point out that the phase-space trajectories of Ref. can be used to find the evolution of the Wigner 

function to the first order in time increments only. 
^ ■ Consider the double Fourier transform of the Wigner function W^(q, p,t): 

^ : W{q, P.r)=J exp(ls(p - p'))T4^'(q, p', r). (1) 

We restrict ourselves to Hamiltonian functions of the form 

H{<l,p) = ^p' + V{<i). (2) 
Using representation ([1]), the quantum Liouville equation 

i> : „ 

—W{q,p,T) = -W{(l,p,T)AH{q,p) (3) 



[ where A is the Moyal bracket 0, H, @ can be represented in the form 



It is tempting to require 



° ^ + ~ ^^^^^ + f ) - ^(q - |)))w^'(q, p', r). (5) 



This is the first-order PDE. It can be solved using the method of characteristics: 

I^'(q, p', r) = l¥'(q - pr, p', 0) exp(^ dr'(F(q - p(r - r') + |) - ^(q - p(t - r') - |)))). (6) 

We observe that Eq.([5]) and its solution Eq.® depend on s and p. 

Given the dependence on s enters into M^'(q, p',t), equation PF(q, p, r) = PF'(q, p,r) does not hold anymore, 
however, Eq.Q is still valid. The appearance of s in W'{q, p', r) does not violate the equivalence of Eqs.Q and ([4]). 
However, given the dependence on p enters W^'(q, p',t), Eq.(I3]) becomes distinct from Eq.Q. 
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There is therefore the obvious shift in the definition of W between Eqs.(l) - (5), since W depends in (5) on a 
parameter. This parameter, p, is not strictly related to the same notation in Eq.(l) on which W depends but W 
does not. The arguments leading to Eqs.fUl) and ([6]) are thereby not founded. 

On the other hand, suppose that Eq.® holds nevertheless. Consider the evolution of one-dimensional oscillator 
with potential V{q) — ^q^. The Wigner function of the ground state has the form 



Applying Eqs. ([T]) and ©, we obtain 



W{q,p)^2exp{-q^/h-p^/n). 



(7) 



W{q,p,T) = 2e^pi-{q + prf /h - {p - qr - pT^ /2)^ /h). 



(8) 



The Wigner function of stationary states does not depend on time. This particular contradiction proves that Eqs. ([5]) 
and ([6]) are incorrect for finite r. 

The expression (|S]) holds however to order 0(r). One can verify, indeed, that Eq.® upon the double Fourier 
transform ffl^ gives the Wigner function correct to order 0(t). 

In Ref. [171 a numerical example is considered to compare the exact evolution of the Wigner function with the 
Wigner function calculated by the infinitesimal step- by-step evolution according to Eq.®. A discrete form of the 
phase-space path integral 13, 2^ is calculated per se. Equation ([6| in its infinitesimal form leads to the representation 
of the path integral discussed by Leaf [12] and Marinov |23| . The first-order solution is therefore useful to obtain the 
correct numerical results [l7j . 

The role of trajectories entering the phase-space path integral is distinct from the role of characteristics in phase 
space. Particle trajectories of the de Broglie - Bohm theory depend on quantum states. In order to find them 
one has to know the wave function. The phase-space trajectories appearing as the Weyl's symbols of the Heisenberg 
operators of canonical coordinates and momenta are currently the only known objects which can be treated as quantum 
characteristics. They comprise the complete information on the evolution of quantum systems and can be used to 
find the evolution of the Wigner function for finite time intervals. 

The authors wish to acknowledge useful discussions with Prof. Amand Faessler. 
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